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SECTION  1 
INTRODUCTION 


Calculations  of  the  source  region  EMP  response  of  very  long 
buried  and  elevated  cables  have  generally  been  performed  by  the  use  of  the 
transmission-line  approximation  to  the  full  electromagnetic  scattering 
problem.  (For  examples,  see  References  1-2.)  That  approximation  is 
useful  when  the  response  is  dominated  by  an  approximately  transverse 
magnetic  and  transverse  electric  wave  that  is  guided  along  the  cable  (a 
TEM  mode).  Usually,  if  a  guiding  structure  supports  a  TEM  mode,  that  mode 
dominates  the  response  of  the  structure  and  the  transmission-line  theory 
applies.  This  appears  to  be  the  case  with  insulated  antennas  and  wires 
since  transmission-line  theory  predicts  their  response  very  accurately. 

The  opposite  situation  exists,  however,  with  a  bare  wire  in  the  soil  or 
air.  Although  such  a  wire  supports  a  guided  wave  that  is  approximately  a 
TEM  mode,  it  is  very  difficult  to  excite.  Indeed,  if  the  wire  were  a 
perfect  conductor,  an  infinite  amount  of  energy  would  be  required  to 
excite  the  TEM  mode. 

The  response  of  a  perfectly  conducting  wire  to  a  delta-function 
voltage  generator  on  the  wire  was  investigated  in  Reference  3.  It  was 
found  that  the  results  could  be  described  by  equations  similar  to 
transmission-line  theory  far  from  the  source  only  if  the  inductance  term 
was  taken  to  be  a  function  of  distance  from  the  source.  The  reason  for 
this  behavior  is  that  the  source  does  not  generate  any  modes  in  the  wire. 
Instead,  it  generates  a  wave  known  as  a  space  wave  that  radiates  out  into 
space.  The  transmission-line  equations  do  not  apply  to  such  waves. 


The  situation  is  much  more  complicated  when  the  wire  has  finite 
resistance.  In  that  case,  the  source  can  excite  both  a  space  wave  and 
electromagnetic  modes.  One  of  those  modes  is  approximately  a  TEM  mode 
that  can  be  treated  by  transmission-line  theory.  The  only  way  to  deter¬ 
mine  the  relative  importance  of  these  various  waves  is  to  solve  the  full 
electromagnetic  scattering  problem.  This  is  the  purpose  of  the  present 
study . 


In  this  report,  the  coupling  to  an  infinitely  long,  finitely 
conducting  wire  in  homogeneous  material  is  examined.  An  exact  integral 
expression  is  obtained  for  the  current  generated  in  the  wire  by  a  voltage¬ 
generating  source  localized  on  the  wire.  The  contributions  of  the  space 
wave  and  the  individual  modes  are  isolated  by  applying  contour-i ntegration 
techniques.  The  contribution  of  the  space  wave  arises  as  a  branch-cut 
integral  and  the  contributions  of  the  various  modes  rise  as  residues  of 
poles.  Roth  the  branch-cut  integral  and  the  residues  are  calculated 
numerically  to  compare  their  relative  importance.  In  addition,  analytic 
approximations  for  the  space-wave  current  are  obtained  for  points  very 
near  and  very  far  from  the  source. 

Numerical  calculations  were  made  for  highly  conducting  lines 
(where  the  results  reduced  to  those  obtained  in  Reference  3)  and  for  Tines 
with  resistances  representative  of  buried  communication  and  power  systems. 
It  was  found  that  the  TEM-type  mode  is  not  excited  in  wires  with  typical 
resistances  at  the  frequencies  of  interest.  Hence,  the  space  wave  is  the 
dominant  contributor  to  the  current  in  these  wires  and  consequently,  stan¬ 
dard  transmission-line  theory  is  not  applicable.  Moreover,  it  was  found 
that,  with  representative  resistances,  the  current  produced  is  quite 
different  from  the  current  that  would  be  produced  on  a  perfectly  conduct¬ 
ing  wi  re. 
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SECTION  2 

INTEGRAL  REPRESENTATION  OF  THE  CURRENT 


The  problem  of  the  current  excited  by  a  delta-function  voltage 
generator  on  the  surface  of  a  straight  wire  has  been  treated  in  the  case 
of  a  perfectly  conducting  wire  by  specifying  the  boundary  value3*4 

E2  =  V06(z) 

on  the  surface  of  the  wire  where  Ez  is  the  component  of  the  electric 
field  in  the  direction  of  the  wire.*  When  the  wire  has  finite  conduc¬ 
tivity,  however,  we  cannot  specify  the  value  of  the  field  anywhere  in 
finite  space.  The  problem  is  best  treated  as  a  scattering  problem  with  a 
specified  field  ET  incident  on  the  wire  generating  a  scattered  field 
Es  that  is  determined  by  the  boundary  conditions  at  the  surface  of  the 
wire.  The  incident  field  is  the  field  radiated  by  the  generator  if  the 
wire  is  not  present.  Our  first  task,  then  is  to  determine  that  field. 

In  a  material  with  finite  or  zero  conductivity,  a  voltage  gener 
ator  produces  a  discontinuity  in  the  electric  field  across  the  generator. 
The  magnetic  field  can  be  continuous.  In  this  study,  the  generator  is 
specified  to  have  these  properties  in  such  a  way  that  the  current  it 
produces  in  the  wire  reduces,  in  the  limit  of  infinite  wire  conductivity, 
to  the  result  already  known  for  the  current  produced  by  a  voltage  genera¬ 
tor  on  a  perfect  conductor.  The  same  generator  has  been  applied  by  Wait5 
and  by  Wait  and  Hill6  to  study  waves  produced  on  cylindrical  rods.  (Thos< 


*  This  problem  was  introduced  into  the  EMP  literature  in  Reference  3.  li 
has  a  long  history  in  the  antenna-theory  literature.  See  Reference  4 
and  the  references  therein. 


studies  differ  from  ours  in  that  they  apply  an  approximate  boundary  condi¬ 
tion  at  the'  rod  surface).  The  generator  can  be  thought  of  as  a  toroid 
with  very  large  current  in  a  wire  very  tightly  wound  around  a  circular 
ring  that  encircles  the  wire.  It  can  be  modeled  mathematically  by  a  ring 
of  magnetic  current.  We  take  the  ring  to  lie  in  the  plane  z  =  U  and  to 
have  radius  b  which  is  larger  that  the  radius  a  of  the  wire.  Later,  we 
set  b  =  a  to  put  the  source  on  the  wire. 

Consider  a  monochromatic  electromagnetic  field  with  angular 
frequency  w  in  an  isotropic  homogeneous  material.  It  can  be  expressed  in 
terms  of  axial  Hertz  vectors  in  the  form7 


E_(r,<{>,z)=VxVx[U(r,<|>,z)lz]  +  ioupV xCV ( r,<f>  ,z  )lz]. 


H(r,<t>,z)=Vx[nU(r,^  ,z)lz]+VxVx[V(r,4>  ,z)lz]  (lb) 


in  SI  units  where  the  scalar  potentials  satisfy  the  inhomogeneous  Helm¬ 
holtz  equation 


(v2+k2 )  U(r  ,4>  ,z)  =  f..(r,<j>  ,z)  , 


(v2+k2  )  V( r »4> , z )  =  f w ( r ,z) 


with  source  functions  fy  and  fy  that  are  nonzero  only  on  sources  or  on 
boundaries  separating  different  materials.  The  time  behavior  e  has 
been  suppressed  in  (1)  and  [2)  and  k  and  n  are  defined  by 


2 

k  =  iyain  , 


(3) 


where  jj ,  e,  a  are  respectively  the  permeability,  permittivity  and 
conductivity  of  the  medium.  We  apply  the  cylindrical  coordinate  system 
illustrated  in  Figure  1. 


m>£2>02  r 


Figure  1.  Magnetic  current  loop  around  a  wire. 

It  follows  from  (1)  that  the  axial  components  of  the  fields  are 

given  by 

E  (r,<|>,z)  =  (k2  +  .?_)  U(r,$,z)  ,  (5) 

2  3z2 

H  (r,4>,z)  =  (k2  +  1_)  V(r,  <t>,z).  (6) 

2  3z 2 

Since  the  generator  produces  an  axial  electric  field  but  no  axial  magneti 
field,  we  assume  that  it  is  a  source  of  the  potential  U  but  not  of  poten¬ 
tial  V.  Since  the  scattering  by  the  wire  does  not  introduce  any  axial 
magnetic  field  either,  we  take  V  to  be  zero  for  both  incident  and 
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scattered  fields.  Moreover,  from  the  symmetry  of  the  problem,  the  poten¬ 
tial  U(r,<}>,z)  is  independent  of  41 .  For  such  fields,  (1)  reduce  to 


: 

Er(r,z)  -»2u<r'z»  , 

3r3z 

(7a) 

r 

E.(r,z)  =  0  , 

4> 

(7b) 

- 

•  ' 

E _(r,z)  =  (k2  ^  )  U(r,z)  , 

3  Z2 

(7c) 

Hp(r,z)  =  0  , 

( 7d ) 

i 

H.(r.z)  -  -n  3U(r’2)  , 

<P  dr 

(7e) 

y 

Hz(r,z)  =  0 

(7f ) 

a 

and  (2a)  reduces  to 

2 

[-  (r| — )  +  d~—  +  k2]  U(r,z )  =  fy(r,z)  . 

r3r3raz2  U 

(8) 

We  now  Fourier  transform  (7)  and  (8)  with  respect  to  z  using  the 
conventions 


A(h)  =  /  A(z)e"ihzdz  , 


A(z)  =  1-  /  A  (h)eihzdh  , 

^TT  wn 


(9a) 


(9b) 
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to  obtain  (for  the  nonzero  components) 


E~r(r,h) 

=  3U(r,h) 

3r 

(10a) 

Ez(r,h) 

=  a2U(r,h)  , 

(10b) 

H.  (r,h) 

,  _n  3U(r,h) 

(10c) 

r 

3r 

[I  —  (r 
r  3r 

+  a2]  U(r,h)  =  f.,(r,h)  , 

3r  u 

(11) 

where* 


We  are  now  ready  to  solve  for  the  field  radiated  by  the 
generator  without  the  presence  of  the  wire.  We  use  a  superscript  i  for 
this  field  to  indicate  that  it  is  the  incident  field  in  the  scattering 
problem.  A  subscript  2  is  used  on  the  material  parameters  e  and  a  and  the 
corresponding  quantities  k , n, a  def i ned  in  (3),  (4),  (12)  to  distinguish 
them  from  the  same  quantities  within  the  wire  which  are  indicated  with 
subscript  1. 

Since  the  source  function  for  the  potential  U1 (r,z)  is  zero 
everywhere  except  on  the  ring  of  radius  b  in  the  plane  z=0,  we  have 

fy (r,z )  =  ^5(r-b)6(z)  ,  (13a) 


*  For  real  h,  the  branch  of  the  square  root  in  (12)  is  chosen  so  that  the 
imaginary  part  of  a  is  positive  or  zero.  The  branch  used  for  complex  h 
is  defined  later. 


and  therefore  have 


f  (r,h)  =  -  6(r-b)  .  (13b) 

U  r 

Hence,  for  r*b,  (11)  is  Bessel ‘s  equation  of  order  zero  with  solutions 
U1  (r,h)  =  CiH^  ^ (a2r)  for  r>b  ,  (14a) 


U7(r,h)  =  C2J0(a2r)  for  r<b  ,  (14b) 

where  the  Hankel  function  in  (14a)  was  chosen  to  insure  outgoing  waves  at 
r="  for  positive  <*2  and  the  Bessel  function  in  (14b)  was  chosen  to  insure 
a  finite  result  at  r=0. 

The  ratio  of  the  constants  in  (14)  can  be  determined  by  requiring 
the  continuity  of  H(j)(r,h)  across  r=b.  According  to  (10c),  this  means  that 
the  r  derivative  of  U(r,h)  must  be  continuous.  Since  the  derivatives  of 
the  Hankel  and  Bessel  functions  with  respect  to  their  argument*  satisfy8’ 

fio^(«r)  =  -H^^(ar)  ,  (15a) 

J  0 ( ar )  =  - J 1 ( ar )  ,  (15b) 

we  find 

C2  H!(1)(a2b) 

—  =  -—-4--  .  (16) 

Ci  Ji  a2b 


*  We  use  a  dot  over  a  function  to  indicate  its  derivative  with  respect  to 
its  argument. 


i9 


-V-V-V- v 
A  -W.' 


5* 


•  ►  »  •  •  «  •  m  *  _ 

-V-Va  ’  -**  .*>  .*■  -*»-*»  .‘•-.vlv . 
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Hence,  (14)  can  be  written 

U1  ( r , h )  =  C0Ji(o2b)  H0 ^ 1 ^ (a2r)  for  r>b  ,  (17a) 

(r,h)  =  C0H1^1^(a2b)  J 0 ( ot2r )  for  r<b.  (17b) 

We  can  solve  for  the  constant  Co  in  terms  of  the  discontinuity 
in  U^r.h)  across  r=b,  but  that  is  unnecessary.  We  determine  it  later  in 
terms  of  the  voltage  V0  across  the  generator.  Also,  we  could  now  assemble 
the  expressions  for  the  incident  electromagnetic  field,  but  again,  there 
is  no  point  in  doing  so.  It  is  most  convenient  to  solve  the  scattering 
problem  in  terms  of  the  potentials  in  Fourier  space. 

To  that  end,  consider  an  infinitely  long  straight  wire  with 
circular  cross  section  of  radius  a.  The  axis  of  the  wire  lies  along  the  z 
axis.  The  permittivity  and  conductivity  are  e!  and  respectively  in  the 
wire  and  £2  and  °2  outside  the  wire.  The  permeability  is  taken  to  be  u  in 
both  materials.  We  take  the  radius  b  of  the  generator  to  be  slightly 
larger  than  a  so  that  its  field  is  incident  on  the  wire  from  the  outside. 
After  solving  the  scattering  problem,  we  can  let  b+a  so  that  the  generator 
is  on  the  surface  of  the  wire.  The  geometry  is  shown  in  Figure  1. 

The  incident  field  excites  a  scattered  field  with  potential 
function  Us(r,z)  that  propagates  outward  away  from  the  wire  in  the  outer 
material.  The  total  field  inside  the  wire  has  potential  function 
Uw(r,z).  Both  Us(r,z)  and  Uw(r,z)  satisfy  the  Helmholtz  equation 
(2a)  with  f|j=0,  and  with  k  replaced  by  ki  in  the  inner  material  and  by 


k2  in  the  outer  material.  Hence,  the  Fourier  transforms  Us(r,h)  and 

Uw(r,h)  of  U$(r,z)  and  Uw(r,z)  respectively  satisfy  Bessel's  equation 
with  a  defined  in  (12)  replaced  by  a!  =  /k2  -  h2  in  the  wire  and 

I  2  2  *  ^ 

a2  ~  \]k2~ h  outside  the  wire.  Hence,  we  must  have 

Us(r,h)  =  AH 0 ^  \a2r)  for  r>a  ,  (18a) 

Uw(r,h)  =  BJ0(a1r)  for  r<a  ,  (18b) 

where  the  particular  solutions  to  Bessel  equations  were  chosen  in  the  same 
way  as  in  (14). 

The  constants  A  and  B  in  (18)  can  be  determined  by  applying  the 
boundary  conditions  that  the  tangential  components  of  the  electric  and 
magnetic  fields  must  be  continuous  at  the  surface  of  the  wire.  According 

to  (10),  this  means  that  a2U(r,h)  and  ryi-i r ih)  must  be  continuous  there, 

3r 

where  U(r,  h)  is  the  Fourier  transform  of  the  total  field.  Outside  the 

"*  ~i  ~c 

wire,  U(r,  h)  is  U  (r,h)  +  U  (r,h)  whereas  inside  the  wire  it  is  U  (r,h). 
Hence,  with  the  use  of  (17b)  and  (18),  the  boundary  conditions  become 

C o «2 ^ 1  ^  )0 o ( qz3 )+Aoi2H [|  )(«2a)  =  Ba1J0(a1a)  ,  (19a) 

Con2H\  ( azb )  *2^  o  ( a2d )  +A  n2q2H  q  (q2a)  =  8njqjJg(qja).  (19b) 

With  the  application  of  (15),  (19)  becomes 


*"Ari 201 2 ^  1  (a2^ )  *  Bn^ai  (a^  a)  ”  (a2  b )Ji  (a2 d )  •  ( 2 Ob ) 


Solution  of  this  pair  of  simultaneous  linear  equations  yields 

n  \  oia2 Jo  (°2  a  )^i  (ai  a )  "“l  Ji  (“2  a ) (ai  a)d2 

A  =  -CqHj  (02b)  -  , 

<*2*1 1  Hq  ^  ^  (a2  a  (a^  a ) -0^02  ^  (<*2  a )  Jq  (a^  a ) 

(21a) 

u  r  a2  u  O)/  u\  ^1  ^  ^  (a2a)Jo  (a2a  )“^0  ^  ^  (<*2a  )^1  (“2 a  ) 

B  =  -C0P2  —  Hi  (a2  b )  - 

al  (1  )  n  \ 

0l2r'  1  Ho  '  (®2  a ) ^1  (al  a )  -<*1^2  Hi  ;  (a2  a )  Jq  (aj  a ) 

(21b) 

The  numerator  in  (21b)  is  a  Wronskian  known  to  satisfy9 

H1(1)(Z)J0(Z)  -  H0(1  )(Z)J1  (Z)  =  J_  .  (22) 

it  i  Z 


Hence,  (21b)  becomes 


B  =  - 


200(12^2 


Hx  (1  )(a2b) 


1  Ql2r|  i  Hq  ^  («2a  )^i  (al  a  )  "alr>2  Hi  ^  ^  (a2  a )  Jq  («i  a ) 


(l  ), 


•  (23) 


2  2 


According  to  (3)  and  (4),  02/01  =  k2/ki*  Hence,  (23)  can  be  rewritten 


B  =  - 


2C 0a2k2 


Hi  (1  }(a2b) 


ai  a2kiH0^  ^  (a2a  )Jj  (ai  a ) -aj  k2Hj  ^  ^  (02  a  )Jg  (04  a  ) 


.  (24) 


We  now  have  everything  we  need  to  calculate  all  field  quanti¬ 
ties.  Since  our  primary  interest  is  in  the  current  generated  in  the  wire 


we  will  concentrate  on  the  field  inside  the  wire.  From  (18b)  and  (24), 
the  Fourier  transform  of  the  potential  of  the  total  field  inside  the  wire 
is 

2Coa2k2  Hi  ^  )(a2b)J0  (cq  r) 

Trial  a  2  (1),  ,  ,  ,  2  (1) 

a2kiHo  (<*2a)Ji  (aia)-a1k2H1  (a2a)J0(a1a) 

(25) 

The  Fourier  transform  of  the  z  component  of  the  current  density  satisfies 

J2(r,h)  =  OiE^r.h)  =  ai01Uw(r,h)  (26) 

and  the  Fourier  transform  of  the  total  current  flowing  along  the  wire  is 
the  integral  of  J(r,h)  over  the  cross  section  of  the  wire 

a 

I(h)=  2n j  J  (r,h)  rdr.  (27) 

o  z 


UW(r,h)  =  - 
for  r  <  a. 


Hence,  by  combining  (25)  -  (27)  and  employing10 


/  rJ0  (a!  r  )dr  =  £-  Ji  (c^r) 
o  1 


"  —  Ji  (aia)  » 


(28) 


we  obtain 


1(h)  =  - 


4C  Q&2  ^2a 


Hi  (  )(a2b)  (at  a  ) 


i  2  (1  )  2  (1  ) 

a2kiHo  (<*2a)Ji  (04  a)-ai  k2Hi  (a2a)Jo(aia) 


The  desired  integral  expression  for  the  total  current  I(z)  flowing  along 
the  wire  follows  immediately  if  we  take  the  inverse  Fourier  transform  of 
(29)  according  to  (9b).  The  result  is 


2*20! 

,<z>  ■  -  —  1  - 2  (T) - 2 — (T) - dh  • 

a2kiHo  (a2a)J  1(a1a)-a1k2H !  (a2a)J0(a1a) 

(30) 

The  expression  in  (30)  is  similar  to  the  formula  applied  by  Wait 
and  Hill6  to  treat  the  current  generated  on  a  metallic  rod  by  a  magnetic 
ring  current.  Their  formula  was  derived5  by  assuming  an  approximate 
impedance-type  boundary  condition  without  specifying  what  the  impedance 
actually  is.  Such  boundary  conditions  are  known  to  be  useful  when  one  of 
the  materials  ha:  a  large  but  finite  conductivity.5  If  we  replace  the 

factor  outside  the  integral  sign  in  (30)  by  so  that  I(Z)  includes 

2 

the  displacement  current  or  if  we  approximate  k2  in  the  first  term  in  the 
denominator  by  iuujc^,  then  (30)  becomes  equivalent  to  the  formula  in  Wait 
and  Hill6  if  we  take  the  series  impedance  of  the  wire  to  be 


l  -  aiJo(aia) _ 

s  2nan1J1(a1a) 

and  our  constant  Cq  to  be 

C0  =  ~  K  ,  (32) 

2ct2 

where  K  is  the  strength  of  the  magnetic  ring  current.  The  expression  for 
the  series  impedance  of  a  wire  given  in  (31)  is  one  that  has  been  used  in 

the  EMP  literature  when  the  definition  of  Zs  is  taken  to  be 

2 

Ez(a,z)/I(z)  for  axi symmetric  TM  fields. 


We  can  express  Cq  in  terms  of  the  voltage  V0  across  the  voltage 

generator  by  relating  (30)  to  the  known  result  for  the  current  produced  on 

a  perfect  conductor  by  a  voltage  generator.  To  that  end,  we  let  b=a  so 

that  the  source  is  on  the  wire,  and  consider  (30)  in  the  limit  as  oi*“. 

In  that  limit,  the  second  term  in  the  denominator  in  (30)  is  negligible 

2 

compared  to  the  first  and  oj  /ki+l/(ipoj ).  Hence,  we  have 

iim  -  2ka  r  C°Hl^  ^a2a)  i  hzHh  i 

lim  I(z)  = - J  - — -  e  dh.  (33) 

ai-  irpw  _  Ho(1)(a2a) 

This  corresponds  to  the  known  expression  in  this  limit3*4  if  we  take 


Comparison  of  (34)  and  (32)  shows  that  the  magnetic  ring  current  of 
strength  K  is  a  voltage  generator  of  strength  V0 . 

The  current  can  now  be  expressed  terms  of  the  /oltage  V0  even 
when  the  generator  is  not  on  the  wire  by  changing  a  to  b  in  (34)  and 
substituting  into  (30) 


I(z)  =  -  V0k2a1b  f 


Hj  (a2b  )Jj  (ai a  )e1  z 
2(1)  2(1) 
a2^i^o  (<*2 a ) Ji  (ex i  3 )  — cx i  kz H i  (a2 a ) Jq  (ai  a ) 


*  m  *  m  “  •  -  •  '  .  '  j'*  ,’*  . 


Equation  (3b)  presents  an  exact  integral  expression  for  the  total  current 
generated  by  the  source  as  a  function  of  distance  along  the  wire.  Unfor¬ 
tunately,  it  is  difficult  to  determine  the  properties  of  the  current  from 
the  integral  in  its  present  form.  In  the  next  section,  we  transform  the 
integral  into  a  form  which  is  more  amenable  to  asymptotic  and  numerical 
analysis.  Moreover,  the  transformed  result  provides  new  insight  into  the 
physics  of  the  current  flow. 
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SECTION  3 


TRANSFORMATION  OF  THE  CURRENT  INTEGRAL 


In  this  section,  we  transform  the  integral  in  (35)  for  positive 
z  by  completing  the  contour  of  integration  in  the  upper  half  plane  and 
applying  Cauchy's  theorem  of  the  residues.*  To  that  end,  we  rewrite  (35) 
in  the  form 


l(z)  =  -V0k2o1bP(z)  , 


(36) 


where 


P(z)  =  /  L (h )dh  , 


(37) 


1(h)  = 


h( 1 )(a2b)J1(ala)eihz 


2  /i  \  2  I  1  \  * 

aik2  H\  ,(a2a)J0(a1a)-a2lc  ^  ’  ( a2a )  J  x  (a1  a  ) 


(38) 


and  investigate  the  analytic  properties  of  the  integrand  L(h)  in  the  upper 
half  of  the  complex  h  plane. 


We  apply  a  procedure  that  is  standard  in  the  theory  of  guided  waves. 
E .g. ,  see  Reference  1 1 . 


20 


3.1 


MULTIVALUED  FUNCTIONS  AND  BRANCH  CUTS 


The  integrand  L(h)  involves  multivalued  functions  04,  a2,  and 
^(Z)  which  must  be  made  single  valued  by  specifying  the  branches  to  be 
used.  For  a! ,  we  choose  the  branch  cut  to  be  the  curve*  in  the  complex  h 
plane  where**  a'i=0.  This  branch  cut  is  shown  as  a  dashed  line  ending 
at  the  branch  points  h=±kj  in  Figure  2.  We  take  04  >0  on  the  top 
Riemann  sheet  and  a"  _<0  on  the  bottom  Riemann  sheet.  For  a2,  we  take  the 

branch  cut***  to  be  the  vertical  lines  terminating  at  the  branch  points 
h  =  ±k2  that  are  shown  in  Figure  2  as  wiggly  lines.  The  dashed  curves 
terminating  at  these  branch  points  are  curves  where  a2=0.  In  the  shaded 
region  between  each  of  these  curves  and  its  closest  branch  cut,  “2  is 
negative  on  the  top  Riemann  sheet.  Elsewhere  in  the  complex  plane,  <*2  is 
positive  on  the  top  Riemann  sheet.  It  follows  from  these  definitions  and 
the  original  definition  of  a  given  in  (12)  for  real  h,  that  the 
integration  contour  in  (37)  is  on  the  top  Riemann  sheets  of  both  04  and 
a2- 

(  1  ) 

For  the  Hankel  function  Fig  (Z),  we  choose  the  branch  cut  to  lie 
along  the  negative  half  of  the  imaginary  axis  in  the  complex  Z  plane. 

This  differs  from  the  standard  definition  of  H0(1)(Z)  which  uses  the 
negative  real  axis.  Our  branch  of  the  Ha nkel  function  is  chosen  so  that 
it  is  equal  to  the  standard  Hankel  function  in  the  upper  half  of  the  Z 
plane.  The  branch  cut  fo,  the  Hankel  function  H0(')(a2a)  occurs  only  on 
the  lower  Riemann  sheet  of  <*2.  On  that  Riemann  sheet,  it  lies  along  the 
dotted  curve  shown  in  Figure  2. 

*  This  branch  cut  is  frequently  called  the  Sommerfield  or  the  fundamental 
cut.  See  Reference  11  for  a  discussion  of  the  characteri sti cs . 

**  We  use  the  notation  that  Z'  is  the  real  part  and  Z"  is  the  imaginary 
part  of  any  complex  quantity  Z. 

***  This  branch  cut  is  called  the  vertical  branch  cut. 


Figure  2.  Conplex  h  plane. 


Even  though  the  point  h=k L  is  a  branch  point  of  alf  it  turns  out 
that  it  is  not  a  branch  point  of  the  integrand  L(h).  To  see  that,  we 
temporarily  write  L(h)  as  L(h,a1)  and  examine  the  value  of  the  integrand 
on  the  opposite  Riemann  sheet  of  at,  by  evaluating  L(h,e1T1a1),  the  result 
obtained  when  the  sign  of  04  is  changed  in  the  right-hand  side  of  (38). 

We  have 


L(h,eniai)  = _ Hi1)(a2b)J1(e”1a1a)elhz _ _ 

e1Ti  “ik 2H i 1  ^(a2a)J0(e1Ti  a],a)  -  a2k  ^ (a2a )J  j (eni  ) 

(39) 

With  the  application  of8 

Jn(eir1z)  =  eniTiJn(z)  ,  (40) 

(39)  becomes 

,  *1  x  _ _ eni  H  [ 1  ^  ( a2b)J  1  ( ot1  a ) _ 

L(h,e  ai)  2  1 1  \  .  2  ni 

eni  04k  2h|  ^(a2a)JQ(a1a)  -  ,(«2a)J1(o1a) 

■  L(h,«i).  (41) 

Hence,  L(h)  has  the  same  value  on  both  Riemann  sheets  of  ai.  This  means 
that  the  point  h=kL  is  not  a  branch  point  of  the  integrand. 

3.2  POLES 

The  integrand  has  poles  where  the  denominator 

0(h)  =  a^k^1  ^(a2a)J0(a1a)  -  «2k  )  ( a2a  )J  l  (ata  )  (42) 


is  zero.  The  equation 
0(h)  =  0 


(43) 


is  a  complicated  transcendental  equation  that  can  not  be  solved 
analyti cal ly . 


For  the  case  of  interest,  however,  it  is  easy  to  obtain  a  rough 
idea  of  the  locations  of  the  poles  in  the  complex  h  plane.  Since  the 
conductivity  of  the  wire  is  much  greater  than  that  of  the  outer  material, 
we  have  |ki|>>  k 2 1  *  Tins  means  that,  in  order  to  have  the  second  term  in 
(42)  be  as  small  in  magnitude  as  the  first  term,  h  must  be  near  a  zero 
of  o2Hg  ^ (a2a)J 1 ( 04a ) .  Since  our  branch  of  Hq  ^(a2a)  does  not  have  any 
zeros*,  the  poles  must  be  located  near  h=k2  (the  zero  of  ot2)  and/or  the 
zeros  of  Ji(«ia).  Since  the  zeros  of  J  i ( Z )  are  all  real,  all  poles  that 
are  not  near  h=l<2,  must  lie  near  the  04  branch  cut  (where  04=0)  shown  as 
the  dashed  curves  terminating  at  ±kt  in  Figure  2.  In  the  limit  as  kt 
tends  to  infinity,  the  poles  near  k2  tend  to  k2  whereas  the  others  tend  to 
i  nf i ni ty . 


For  very  large  |h| ,  ki  can 
expression  (12)  for  04  leading  to 

04  =  ih  for  h‘>0  , 

04  =  -ih  for  h 1  <0  . 

In  this  region,  the  large-argument 

jl(Z)  :JL costz  -|I) 

can  be  applied  to  obtain  the  locat 
shows  that  any  poles  present  in  th 
poi nts 

h"  =  ±  (N  +  I)  - 
4  a 

on  the  imaginary  axis  for  large  integral 

*  See  Figure  9.4  of  Reference  8.  The  z 
on  a  different  Riemann  sheet  when  our 
calculations  have  shown  that  no  zeros 
branch . 


2 

be  neglected  compared  to  h  in  the 

(44a) 

(44b) 


approximati on 


ons  of  the  zeros  of  JjJaja).  This 
s  region,  must  be  located  near  the 

(46) 

N. 

ro  shown  in  the  third  quadrant  is 
branch  cut  is  used.  Numerical 
exist  in  that  quadrant  for  our 
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3.3 


COMPLETION  OF  THE  CONTOUR 


For  positive  z,  we  consider  the  integral 

P  (z)  =  /  L(h)dh  (47) 

c 

with  the  same  integrand  as  the  integral  P(z)  defined  in  (37)  but  evaluated 
over  the  closed  contour  c  shown  in  Figure  2.  The  contour  passes  from 
h'=-R  to  h'=R  along  the  real  axis  and  then  along  a  circular  arc  of  radius 
R  into  the  upper  half  plane  until  it  reaches  the  a2  branch  cut.  It  then 
passes  down  along  the  branch  cut,  around  the  branch  point  at  k2 ,  and  back 
up  to  the  same  height  on  the  branch  cut  as  where  it  started.  It  then 
continues  towards  the  left  along  a  straight  horizontal  line  until  h' 
reaches  -k2.  From  there,  it  follows  a  circular  arc  of  radius  R  until  it 
returns  to  the  real  axis.  As  R  tends  to  infinity,  the  real  axis  portion 
of  Pc (z )  becomes  equal  to  our  original  integral  P(z). 


We  begin  by  examining  the  value  of  the  integral  over  the  portion 

of  c  that  is  off  of  the  real  axis.  On  that  portion  of  c,  IhU®  as  R+« . 

2  2  2  11 
Hence,  we  can  neglect  k1  and  k2  compared  to  h  in  the  definitions  of  a! 

and  a2  so  that  satisfies  (44)  and  a2=  -ih  to  the  left  of  the  branch  cut 

and  a2=ih  to  the  right  of  the  branch  cut.  Since  it  was  shown  in  Section 

2.1  that  the  value  of  L(h)  is  unaffected  by  a  change  is  sign  of  04  ,  we 

takea^ih  for  both  positive  and  negative  h‘.  For  large  |h|,  we  can  use 

the  large  argument  approximations8. 


J»(Z)-'^Js1n(Z  -i)  , 


(48b) 


(!) 

Ho  (Z) 


(Z  -f) 

4 


H(x)/7.  -  .  rr  i(z  -  j) 

H1  (Z)  =  -i  y  —  e  4  , 


(48c) 


(48d) 


valid  for  large  jzj  to  approximate  the  functions  in  the  integrand  given  in 
(38).  The  result  is 


L (h >  -  ^eIh(b-a)+ih2K(h), 


(49) 


where 


K(h) 


sin(iha  - 


?k  i  si  n  (iha  -  £)  -  ik^  cos(iha 


(50) 


with  the  top  sign  applying  to  the  right  of  the  vertical  branch  cut  and  the 
lower  sign  applying  to  the  left. 

The  factor  K(h)  is  infinite  at  the  locations  of  the  poles  which 
are  near  points  on  the  imaginary  axis  with  h"  given  by  (46).  In  order  to 
make  sure  that  c  does  not  pass  through  a  pole,  we  choose  R  to  given  by 


R  =  (2N  +  -)  - 
4  a 


(51) 


for  large  integral  N.  We  then  let  R-*-00  in  discrete  steps  by  letting  N+°°. 


■x.v.v.Y- i'.v.; 


■  .  — .  . . 
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Now,  consider  the  portion  of  the  contour  that  passes  from 
h  =  k2  to  h  =  -k2  along  a  horizontal  straight  line.  For  large  enough  N, 
h"  is  approximately  equal  to  R  on  this  line.  Hence, 


h'  =  (2N  +  1)  2L 
4  a 


on  the  contour.  Since  the  poles  are  near  points  that  satisfy  (46),  the 
value  of  K(h)  is  finite  along  the  contour.  Moreover,  since  the  functions 
in  K(h)  are  periodic  functions  of  h"  with  period  of  2it  for  fixed  h',  K(h) 
is  independent  of  N.  Hence,  from  (49),  we  have 


I,  / h \  |  <  M  -h‘‘z 

lL(h)l  "ThT  e 


where  M  is  the  maximum  value  of  ^  j K (h ) | eh  ^b_a^along  the  contour.  Hence, 

the  value  of  the  integral  along  this  part  of  the  contour  tends  to  zero 
exponentially  as  N+“  for  positive  z. 

Next,  consider  the  curved  portion  of  the  contour.  Here,  we  have 
jh  1 1  +0°  as  R+00.  Hence 

cos (i ha  -  tt/4)  =  I  (e-ha-iir/4  +  eha+iTt/4)  =  i  e±(ha+ilr/4)  ,  (54a) 


sin(iha  -  tt/4)  =  i_  (e-ha-iir/4  _  eha+i*/4)  I+l  e±(ha+i  ^/4) ,  (54b) 
2i  2i 


for  large  R  where  the  top  sign  is  used  for  the  right-hand  part  of  the 
contour  and  lower  sign  is  used  for  the  left-hand  part.  Hence, 


/.  *.  *.  \  ■.  *  ,  - 
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K(h)  =  *  -2— 2-  (55) 

k !  +k2 

for  large  R.  Combination  of  (5b)  and  (49)  yields 

L(h)=  +  - 1 — —  exp[  ( - 1  h  '  |  +  ih")(b-a)]e1hz  (5b) 

bh(k1+k2) 

for  large  R.  Since  b-a_>t),  this  is  of  the  form 

L(h)=f (h)elhz  ( 57 ) 

where  |f(h)|  +  U  as  jhl-*-00.  Hence,  it  follows  from  the  usual  proof12  of 
Jordan's  lemma,  that  for  positive  z,  the  integral  over  this  part  of  the 
contour  tends  to  zero  as  R-*». 

Since  the  integral  along  the  portion  of  c  that  tends  to  infinity 
in  the  upper  half  plane  has  been  shown  to  tend  to  zero  as  R+°°,  the  only 
portions  that  are  left  to  contribute  are  the  real  axis  and  the  branch 
cut.  We  have 

Pc(z)  =  P(z)+B(z)  ,  (58) 

where  P(z)  is  the  integral  of  interest  given  in  (37)  and  B(z)  is  the 
integral  around  the  branch  cut.  According  to  Cauchy's  theorem  of  the 
residues12 


Pr(z)=2ni l  R  (z) 

^  i  J 


(59) 
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where  Rj(z)  is  the  residue  of  the  jt*1  pole  and  the  sum  is  over  all  of 
the  poles  in  the  upper-half  plane  on  the  top  Riemann  sheet  of  a2.  Henc 
the  integral  of  interest  has  been  transformed  to  the  form 

P(z)  =  2*il  Mz)-B(z)  (60) 

j  J 

for  positive  z.  It  remains  to  evaluate  the  two  terms  on  the  right-hand 
side  of  (60)  and  to  interpret  them  physically. 

3.4  COMPLEX  MODES 

The  integrand  in  Pc(z)  can  be  written  in  the  form 

L  ( h )  =  ( a2b  )J  i  ( aia  )e / 5  ^  \ 

0(h) 

with  D(h)  defined  by  (42).  Since  the  poles  hj  occur  at  the  zeros  of 
D(h)  so  that 

D(h  )=0  ,  (62) 

the  residue  of  the  pole  is  given  by13 
..( l) ,  .  >  ,  ,  .  ihz 

R  .  ( z )  =  [ill — - ]  .  (63) 

^  D(h)  h=hj 

Since  its  z  dependence  is  of  the  form  eih3z  ,  we  see  that  R  j ( z ) 
represents  a  wave  field  with  axial  wave  number  hj. 


Hence,  (6b)  becomes 


D(h)  »  -h{kHi'1(o,2a)J„U1a)+il|ii[H0(I,(«2a)  -  j0(„ia) 

2 

-k2  aH^  ^ (a2a  )J} (a^a )  ^  (“2  a  )J^  {a\ a )  + 

2 

+kf aHf 1  ^(<x2a)J,  (a1a)  -  H^1  ^(a2a)[  J0  (c^a)  -  ^a^a- -] }  .(67) 

al  U1  ° 


Combination  of  terms  and  application  of 


2  2 
a2  "al 


> 


(68a) 


It  follows  from  the  equation  (43)  for  locations  of  the  poles 
that  the  first  and  last  terms  in  the  square  bracket  in  (69)  are  equal  when 

D(h)  is  evaluated  at  a  pole  hj.  Hence, 


0(n  )  = 

J 


_  "J 


l*v 


4 


ala2 


(ata)  - 


2  (1  ) 

ah  at ^  Hq  (a2  a )  Jq  (a2  a , 
j 


~aot|a2^1  ^  (a2  3  )  <J  i  (a  i  3  )  ]  .  (70) 

where  h  in  al  and  a2  is  to  be  evaluated  at  hj .  Equation  (70)  is  the 

final  exact  expression  for  f)  (hj)  to  be  substituted  into  the  expression 
(63)  for  the  residue  of  the  pole  hj. 


The  properties  of  the  modes  of  a  wire  in  air  have  been  studied 
by  Sommerfeld  (see  Reference  14  for  an  account  of  Somnerfeld's  results). 

He  found  that  there  is  one  mode  with  wave  number  h0  near  k2  on  the  upper 
Riemann  sheet.  He  called  this  mode  the  principal  mode  and  the  others  he 
called  secondary  modes.  The  secondary  modes  are  attenuated  much  more  than 
the  principal  one.  In  fact,  since  the  wavenumber  of  the  secondary  modes 
lie  near  the  dashed  line  in  Figure  2  that  goes  above  lq  in  the  complex 

plane,  h">kj‘  for  all  of  these  modes.  Hence,  they  are  attenuated 

j 

even  faster  than  is  a  plane  wave  propagating  in  the  metal.  This  shows 
that  the  secondary  modes  are  damped  immediately  and  can  be  neglected  after 
a  very  short  distance  along  the  wire.  For  this  reason,  we  now  concentrate 
on  the  principal  mode,  which  we  label  by  j=0. 


Since  h0=k2  and  since  ki  >>  k2 |  ,  the  second  term  in  the  square 
brackets  in  (7(J)  is  negligible  compared  to  the  first  term.  Hence,  the 
expression  (63)  for  the  residue  of  the  pole  becomes 


Ro(z)  =  - 


aia2^i  ^(a2b)elh°z 


ho  (k  r 


2 

k2: 


[  2k 


V1), 

iMo  1 


a2a)-aa1oi2H|  \a2a[ 


(71) 


«ut.i 


V  ‘-i.W 


where  h  in  04  and  «2  is  to  be  evaluated  at  h0.  Next,  we  apply  the 
2  2-2 

approximations  k  1  -k 2  =  and  04  =  to  obtain 


Ro(z)  = 


a2H^  ^(a2b)e''boZ 


h  0k  1  [ 2H  0  1  ^  (  u2a  )  ~A  a2H  1  4 ;  (  “29  )  1 


(*), 


(72) 


Finally,  since  |a2b|  is  very  small  for  the  principal  mode,  we  apply  the 
smal 1 -argument  approximations  of  the  Hankel  functions* 


.P>, 


(73a) 

ir  2i 

-  .  2i 
irZ  ’ 

(73b) 

valid  for  small  j Z j  where  F=eY  and  y=0. 577216  is  Euler's  constant.  The 
result  is 


R0(z)  =  -  - - - -^°Z- - .  (74) 

bh0ki[21n(-^i)  +  1] 

where  a2  is  to  be  evaluated  at  h=h0. 


*  Reference  8,  page  360.  The  factor  1/i  in  the  argument  of  the  logarithm 
in  (73a)  accounts  for  the  contribution  of  Jq(z)  =  1  for  small  Izl.  We 
apply  the  principal  branch  of  ln(z)  which  has  branch  cut  lying  along 
the  negative  real  axis.  With  this  choice,  the  right-hand  side  of  (73a) 
has  the  same  branch  cut  as  the  left-hand  side. 


According  to  (60),  the  contribution  of  R0(z)  to  the  integral 
P(z)  is  2 it i R o  (z ) .  Hence,  it  follows  from  (36)  that  the  contribution 
Ip(z)  of  the  principal  mode  to  the  current  in  the  wire  is 

Ip  *  -  V0k2  oib2TTiR0(z)  (75) 

Application  of  our  approximate  expression  (74)  for  the  residue  Rq(z) 
therefore  gives 

,  ,  2»iVo4oielh<|Z  (76) 

P  h0k^[21n  (Ifia>l] 

2i 

2  _  1 

Since  oj/k  ^(i  uw)“  as  o^®  [see  (3)  and  (4)],  we  see  that  Ip+0  as 
This  means  that  the  principal  mode  is  not  excited  by  the  generator  if  the 
wire  is  perfectly  conducting,*  It  should  be  noted,  however,  that  Ip 
tends  to  zero  rather  slowly  as  tends  to  infinity  because  1  n ( ot2a )  tends 
to  infinity  slowly.  This  means  that  the  contribution  of  the  principal 
mode  to  the  current  might  be  significant  even  with  highly  conducting 
wires.  To  answer  the  question  definitively,  it  is  necessary  to  solve  for 
h0  numerically  in  order  to  obtain  an  accurate  value  of  a2. 

Another  important  question  that  must  be  answered  numerically  is 
whether  the  wavenumber  ho  lies  on  the  upper  or  lower  Riemann  sheet.  For 
the  case  treated  by  Sommerfeld  (see  Reference  14),  it  is  known  that  h0 
lies  on  the  upper  sheet.  As  the  parameters  vary,  however,  h0  can  cross 
the  branch  cut  to  the  lower  sheet.  In  that  case,  the  residue  of  the  pole 
does  not  contribute  to  the  current. 


*  It  can  be  shown  easily  that  the  same  is  true  for  all  of  the  other  modes 
as  well.  This  is  the  reason  that  modes  do  not  appear  in  the  analysis 
in  Reference  3. 


3.5 


SPACE  WAVE 


In  order  to  obtain  a  physical  interpretation  of  the  branch  cut 
integral  B(z),  we  examine  its  asymptotic  behavior  for  large  z.  Since  the 
imaginary  part  of  h  increases  as  we  move  up  either  side  of  the  branch  cut 
from  the  branch  point  k2»  the  integrand  L(z)  decreases  exponentially. 
Hence,  for  large  enough  z,  only  the  vicinity  of  the  branch  point  contri¬ 
butes  to  the  integral.  Near  the  branch  point,  it  is  useful  to  write  D(h) 
in  the  form 


D(h)  =  a1k2H(11)(a2a)J0(a1a)  [l  -  ] .  (77) 

“iMi  (a2a)J0(aia) 

Since  a2=0  at  the  branch  point,  the  second  term  in  the  bracket  is  much 
smaller  than  1.  Hence,  we  can  expand  D_1(h)  in  the  form 

J_:  [aIk22Hj1)(a2a)J0(aIa)]~1[H  +  ...].  (78) 

04k 2H )  '(a2a)J0(aia) 

Neglecting  the  higher  order  terms  and  substituting  into  the  definition 
(38)  of  the  integrand  L(h),  we  find 

L(h)  =  Hi  \ ( a2b)J  1  ( aid )  []+  a2kiHl)  ^(a2a)Ji(aia)  |  e i hz  (79) 

aik2H( 1 ^(a2a)J0(a1a)  ajk^H^1 ^(a2a)J0(a1a) 

Since  |a2|b  and  j«2|a  are  small  near  the  branch  point,  we  apply  the 
smal 1 -argument  approximations  (73)  for  the  Hankel  functions  in  (79). 


The  result  is 


L(h)  =  aJ-t  (°Lij -  [1 

bet  i  ^  2  o  ( aj  a ) 


ak  ^a2  J  i  (a ^  a  ) 

~~2 - 

^2al  ^0  (al  3  ) 


In  I  E5ii)]  eihz 
2i  n 


(80) 


We  next  determine  the  discontinuity  of  L(h)  across  the  branch 
cut  which  results  from  a  change  in  sign  of  a2 .  We  have 


a2 


(81) 


The  first  square  root  in  (81)  is  the  same  on  opposite  sides  of  the  branch 

cut  because  k2+h  does  not  change  as  h  encircles  the  branch  point.  The 

second  square  root,  however,  is  different.  The  argument  of  k2-h  is  -n/2 
on  the  left  side  and  3rr/2  on  the  right.  Hence,  the  argument  of  yj  k2  -h  is 

—tt  /  4  on  the  left  and  3ir/4  on  the  right.  Since  1  nz  =  l  n  |  z  j  +i  (argument  of  z) 

with  -ir  <(argument  of  z)<ir,  the  discontinuity  of  the  logarithm  in  (80)  can 
be  written 


In  ( 


Ta2  a  -j 
2i  J 


-ln(^li) 

left  side 


=  -it  i  . 

right  side 


(02) 


The  corresponding  discontinuity  in  L(h)  from  left  side  to  right  side  is 
therefore 


it i  |-ak1J1  (ega^ j 2 c i h z 
^  k2  Jq  (04  a  )aj 


(83) 


As  a  result,  the  branch  cut  integral  is 


The  integral  B0  can  be  evaluated  exactly.  Let 


h  =  k2+iK  • 
Then  B0  becomes 


B0  =  /  [k2  -  (k2+2i<k2 -<2  )]  e^k2z  e-lCZick 
0 


=  ie1  k2 +  2i k2— )  r  e-<zd< 
8  z  3z  0 


=  ie^k2z(^—T  +  2i k2? — )  A 

77  3Z;  z 


=  ^  elk2z[ 1 -i k 2 z )  . 


Substitution  of  (88)  into  (84)  yives  the  asymptotic  expression 

B(z)  =  —  [ ak  iJ  l  ( ^ia  >  ]  e i k  ?_ z 

b  k?J„(i,a)i,  z3 


valid  for  large  z. 


For  high  frequency  applications,  the  factor  in  the  brackets  in 
(89)  can  be  simplified  by  using  the  approximation  aj  =  kj  and  noting  that 
| k 1 1  a  >> 1  for  typical  wires.  14Hence,  we  can  use  the  large-argument 
approximations  (48a-b)  of  the  Bessel  functions.  Since  k i a  has  large 
positive  imaginary  part,  (48a-b)  can  be  written 

Jo(“ia)  =  - — j; —  e-*  ia-ir/4 )  t  (90a) 

V2^k7a' 


J^a)  =  —  i -  e-1'  (kia-ir/4)  #  (90b ) 

V2^a 


Consequently,  we  have 


J  l  ( aia )  I  .j 

Jq (®ia ) 


(91) 


The  approximate  expression  for  B(z)  then  becomes 


B (z)  :  ^  eik*z 

b  k2  z 


for  large  z  and  high  frequency. 


(92) 


The  contribution  Ir(z)  of  the  branch-cut  integral  to  the 
current  in  the  wire  is  found  by  combining  (89),  (60)  and  (36)  to  obtain 


I 

B 


(2)  «  - 


2jrV0a1 


j-ak  i  ( aja )  j2  l-ik2z  gik2z 

Vo(V>  z" 


for  large  z. 


(93) 


On 

Jv  ) 


The  quantity  Bq  given  in  (88)  is  proportional  to  the  z  component 
of  the  electric  field  radiated  by  a  point  electric  dipole15  located  at  the 
origin  in  a  homogeneous  medium  with  wavenumber  k2-  We  are  led,  therefore, 
to  the  following  physical  interpretation.  The  branch-cut  integral  yields 
a  current  driven  by  the  component  along  the  wire  of  the  electric  field  of 
an  electromagnetic  wave  that  is  radiated  out  by  the  generator.  Far  from 
the  generator,  this  wave  has  the  same  spatial  dependence  as  it  would  have 
if  it  were  being  radiated  in  the  outer  material  without  the  wire  present. 
The  primary*  effect  of  the  wire  is  to  weakly  influence  the  magnitude  of 
the  wave  through  the  multiplicative  factor  in  (84).  Note  that  this  is 
true  even  inside  the  wire,  where  the  wavenumber  is  very  different  from 
k2-  Because  of  these  properties,  the  electromagnetic  field  corresponding 
to  the  branch-cut  integral  is  frequently  called  a  "space  wave."11  It  is  a 
field  that  is  radiated  into  space  rather  than  being  guided  by  the  wire  (as 
are  the  modes).  The  geometrical  factor  (l-ik2z)/z3  occurring  in  (93)  is 
due  to  the  "radiation  spreading"  of  the  wave.  The  reason  this  fall  off  is 
faster  than  the  1/z  behavior  of  a  radiation  field  is  that  the  current  is 
produced  by  the  z  component  of  the  electric  field  which  is  the  radial 
component.  The  radiation  field  that  falls  off  as  1/z  is  the  transverse 
electromagnetic  field  which  does  not  generate  currents  along  the  wire. 

Although  the  approximate  expression  (92)  for  B(z)  is  independent 
of  wire  conductivity,  the  range  of  validity  of  (92)  is  highly  dependent  on 
oj.  This  is  evident  from  the  derivation  of  the  result.  In  order  for  us 
to  be  able  to  apply  the  expansion  used  in  (78),  it  is  necessary  for  the 
second  term  in  the  square  brackets  in  (77)  to  be  much  smaller  than  one  in 
magnitude  over  the  portion  of  the  branch  cut  where  the  integrand  is  not 
negligible.  But  the  magnitude  of  the  second  term  grows  rapidly  with 
increasing  distance  from  the  branch  point  because  | k ! | 2/ j k 2 1 2  is  large. 

*  There  is  another  "effect"  of  the  wire  discussed  in  the  next  paragraph. 


Hence,  z  must  be  large  enough  so  that  the  factor  e^z  causes  the  inte¬ 
grand  to  be  negligible  all  along  the  cut  except  very  near  the  branch 
point.  The  larger  is  the  conductivity  of  the  wire,  the  faster  the  second 
term  in  (77)  grows  and  therefore  the  larger  z  must  be  in  order  for  the 
result  (92)  to  be  valid. 

We  now  obtain  a  rough  estimate  of  the  range  of  validity  of 
(92).  The  second  term  in  (77)  is  nearly  equal  to  1  at  the  locations  of 
poles.  The  pole  h0  is  located  very  near  and  slightly  above  the  branch 
point.*  Hence,  the  second  term  in  (77)  becomes  comparable  to  1  when  the 

point  of  integration  on  the  branch  cut  reaches  the  height  hg  of  the  pole. 

Therefore,  in  order  for  our  expansion  of  (77)  in  the  form  of  (78)  to  be 

,  »  ,  " 
valid,  the  factor  e'noz  must  be  negligible  compared  to  e~K2z.  The  first 

quantity  will  be  down  by  a  factor  e"  compared  to  the  second  if 

( -h g+k •? )z  =  -1.  Hence,  the  asymptotic  result  is  valid  if  z  is  much  larger 

than  the  distance 


z 

c 


h0-k2 


(94) 


This  distance  can  be  very  large  because  h0  can  be  very  close  to  k2.  Since 
ho  approaches  k2  as  the  wire  conductivity  increases,  the  asymptotic  result 
is  valid  only  at  very  large  distances  for  a  highly  conducting  wire.  As 
a t *«=,  the  result  is  not  valid  anywhere.** 


It  is  helpful  to  consider  a  specific  example.  Sommerfeld  has 
shown14  that 


*  This  is  known  from  previous  numerical  studies  of  the  principal  mode. 
See  Reference  14. 


In  this  limit,  the  results  of  Reference  3  and  4  are  valid. 


(95) 


h0  =  k2[l+(6.0+6.4i )10-5] 

for  the  case  of  a  1mm  radius  copper  wire  in  air  with  frequency  of  109  Hz. 
Since  k2  =  0  for  air,  (94)  and  (95)  together  require  that 

i  4 

k2zc  =  1.6x10 


Hence,  z  must  be  much  more  than  2.5xl03  wavelengths  in  air  or  about  720m 
for  (92)  to  be  valid  in  this  case.  Other  numerical  examples  are  presented 
in  Section  3. 

3.6  NEAR  ZONE 

Since  z  must  be  so  large  for  the  asymptotic  result  to  be  valid, 
we  obtain  in  this  section  a  rough  idea  of  the  behavior  of  the  current  much 
closer  to  the  source.  We  wait  until  later  to  specify  how  small  z  must  be 
for  the  approximations  to  be  valid,  but  it  can  be  stated  at  this  point 
that  z  must  be  much  less  than  zc.  On  the  other  hand  we  require  that  z 
be  large  enough  so  that  the  smal 1 -argument  approximations  can  be  used  for 
the  Hankel  functions  over  the  portion  of  the  branch  cut  where  the 
integrand  is  not  negligible. 

On  the  branch  cut,  we  have  h=k2+i<  with  <>0  according  to  (87). 

Since 


41 


(96) 


2 

“2  5  (k2+h)(k2-h)  =  -2ik2< 

for  small  *,  we  see  that  the  arguments  of  the  Hankel  functions  are  small 
if  < << ( 2 j k 2 j b2 ) - 1 .  We  want  the  integrand  to  be  negligible  for  <  outside 
this  range.  To  that  end,  we  require  that  z»zd  where  zd  is  defined  by 

2 

zd  =  2  J  k  2 1 b  .  (97) 

Then,  the  integrand  is  negligible  for  <  too  large  for  the  small  argument 
approximations  to  be  used  because 

jpi  (k2+i  <)zd/e1*<2Zdj  <e-1  (98) 

2  l 

for  k  _>  (2 1 k 2 1 b  )"  .  Hence,  the  region  of  interest  satisfies 

Zd«z<<zc  (99) 


with  zd  and  zc  defined  by  (97)  and  (94)  respectively.  Since 


1  »2b2|k2|  (1001 

h0-k2 


for  highly  conducting  wires  with  small  diameters  and  with  the  source  near 
the  wire,  (99)  covers  a  wide  range  of  distances. 

We  now  employ  the  smal 1 -argument  approximations  (73)  for  the 
Hankel  functions  in  the  expression  (38)  for  the  integrand 


L(h)  = 


aJ  i(aia)e 


i  hz 


bc»|k 2J  g  ( ai  ^  ) "*"3 ba2k  i<J  j  ( aia  )  In  ( — — ) 
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I  J  ,  a  «\  «.  ■  o  o. 


-  '  ' 


(101) 


This  expression  can  be  simplified  further  by  noting  that  when  h  is  near 
k  2 ,  «i=ki  for  | k  i  J  >> | k  2  J -  Then  (101)  becomes 


L(h)  = 


aJ !  (k  iaje1^2 

bk  i  [k  2<J  i(k  ia )+ak  x otj J  i  (k  ja  )ln  )  ] 


The  discontinuity  in  the  logarithm  in  (102)  across  the  branch 
cut  has  already  been  discussed  and  is  presented  in  (82).  It  follows  from 
that  discussion  that  the  argument  of  [ (k2-h )/i  ] 1//2  is  given  by 


arg  \'i  (h-k2)  =  ±  - 


(103) 


on  the  branch  cut  where  the  top  sign  applies  to  the  right  side  and  the 
bottom  sign  applies  to  the  left.  Hence,  the  logarithm  can  be  written 

=  1n(r-a^~4k2+-)  +  ln(VT(h^k7) 

2i  2 

=  in  (££^J^-2)  +  1  n t  1  i 
2  2 


=  In  k2<)  ±1  i  , 


(104) 


where  the  approximation  k2+h  =  2k2  has  been  used.  To  eliminate  the  square 
root  in  (104),  it  is  convenient  to  write 


In (-—•?— )  =  |  [Tn  (A<) )  ±  Hi  ]  , 

2i  2 


(105) 


where 
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where  L  +  (h)  is  the  integrand  on  the  right  side  of  the  branch  cut  and  L_(h) 
is  the  integrand  on  the  left  side  of  the  branch  cut. 

Since  the  contour  of  the  branch-cut  integral  passes  down  the 
right  side  and  up  the  left  side  of  the  branch  cut,  the  integral  can  be 
wri tten 


B(z)  =  i  /  [ I — (h )  -  L+(h)]d<  (108) 

o 


with  h  =  k2+i<.  If  we  introduce  the  constant 


W  =  i  J° (kla) 
Jldqa)  ’ 


(109) 


the  integrand  in  (108)  becomes 


L  -  (h  )  -  U(h)  = 


2ak1iofe'' 


bkl^2  [  -i Wk2  -i  ak^  1  n(A<  )]*  -( ak^K  )' 


(110) 


Hence,  the  branch-cut  integral  is 


2 

B (z )  =  -  2iia-  eik2z  /  _ _ ?  .  (Ill) 

^k 2  0  [Wk 2+ak i <1 n (Ak)  ]  +(ak Lw<) 

We  note  that  for  high-frequency  applications  when  the 
large-argument  approximations  can  be  used  for  the  Bessel  functions,  then  W 
reduces  to  1  according  to  (91).  Otherwise,  W  is  a  constant  that  depends 
on  the  frequency  and  the  properties  of  the  wire. 

For  <  very  near  zero,  the  denominator  in  (111)  is  dominated  by 
the  first  term  Wk2  in  the  square  brackets.  The  second  term  grows  very 
rapidly  with  increasing  k,  however,  because  | k  x | >  > | k  2 1 ,  and  eventually  it 
dominates  over  the  first.  Let  <m  be  the  smallest  value  of  <  such  that 
Wk2  can  be  neglected  in  (111)  for  Then  B(z)  can  be  approximated  by 


B  (z)  =  Bm(z)  +  Bp  (z ) 


where 


(112) 


B  (z)  =  -  eik2z  j<m 

m  bk2 


-<z 


<e 


d< 


rr; 


[Wk2+ak  x<ln (A<) ]  +  (ak  x<n) 


(113) 


B  (z)  =  -  ^T- 
p  bktk2 


ik  jZ 


I 

'Sn 


e_KZd< 

<  ( 1  n  2  (A  < )  +  if2"  ] 


(114) 


It  is  difficult  to  learn  much  about  Bm(z)  without  employing 
numerical  techniques.  Two  points  can  be  made,  however.  First  we  note 
that  for*  zq<<z<<l/Km,  the  z  dependence  of  Bm(z)  arises  only  from 
the  factor  eik2z.  Hence,  Bm(z)  does  not  decrease  due  to  the 
geometrical  spreading  of  the  space  wave.  Secondly,  we  note  that  since 
<m+0  as  °i>00,  Bm(z)  tends  to  zero  as  the  wire  tends  toward  a  perfect 
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Since  <m>ho-k2,  this  restriction  requires  that  z<<zc. 
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conductor.  This  means  that  Bm(z)  may  be  small  compared  to  Bp(z)  for 
highly  conducting  wire. 

We  now  turn  our  attention  to  Bp(z).  The  integral  in  (114)  is 
the  same  as  occurs  in  the  treatment  of  a  perfect  conductor  in  Section  4  of 
Reference  3  except  that  (a)  our  factor  A  defined*  in  (106)  is  complex 
whereas  it  is  purely  imaginary  in  Reference  3  and  (b)  the  lower  limit  of 
our  integral  is  <m  whereas  it  is  zero  in  Reference  3.  As  a  result,  we 
can  approximate  Bp(z)  by  the  procedure  used  in  Reference  3  with  minor 
modifications. 

To  that  end,  consider  the  integral 


I(z)  •  I 


-<Z. 

e  dK 


<m  f[lnZ(A<)+Tt" 


(115) 


Because  of  the  factor  e~<z,  the  integrand  is  negligible  for  <  larger 
than  some  number  a  which  is  small  enough  that  | In (A*) | >>t  for  <<a.  Hence, 
the  term  in  the  square  brackets  can  be  expanded  in  a  binomial  series  to 
obtain 


00 

I(z>  %Si 


(116) 


where 

I  =  (-i)"^2"-2  f  e  —  . 

n  <ln<-n(A<) 


(117) 


* 


Note  also  that  our  A  is  defined  differently  than  the  A  of  Reference  3. 


I 


We  now  change  the  integration  variable  to 
C  =  -ln(A<) 


(118) 


to  obtain 


exp [-  L  e~^ ] 

I  =  (-l)n+lir2n-2  / _ A__  dc 


(119) 


where  is  a  horizontal  straight  line  a  distance  arg(A)  below  the  real 
axis  in  the  complex  c  plane  starting  at  -°°-iarg  A  and  ending  at 
-ln(A<m). 

The  exponent  -  —  e~?  is  real  and  negative  all  along  c^.  Hence, 
A 

the  integral  over  the  portion  of  the  contour  to  the  left  of 


Co  =  -In (A/z ) 


(120) 


can  be  approximated  by 


C0  exp [-  f  e‘c  ] 

l(-)=  (-l)n+lw2n-2  /  - A_ - dc 

n  -°°-iargA  C02n 


(121) 


because  the  integrand  decays  so  rapidly  as  C  moves  a way  from  C0  along  c^ 
to  the  left. 
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In  the  integral  to  the  right  of  Sq,  the  exponential  term 

exp[-  —  e"^]  is  approximately  equal  to  one  except  in  the  immediate 
A 

vicinity  of  Cq.  Hence,  we  approximate  this  integral  by 


exp [-  —  e  ^]-l 

(-1  )n+lTr2n-2  j-ln(A<m)  j-j _  +  _ A _ ] 

Ao  42n  4o2n 


=  (-l)n+17r2n"2{-  J_L_1_  ]2n_1  +  J _ 1 


2n-l  -1 n (A ^ ) 


2n-l  £Q2n-l 


!  In  \  -ln(A<  ) 

ln(AK )  ,  ,  m  _ 

+ - 2-  +  — - -  +  — —  /  exp[-  -  e“c  d?}.  (122) 

Co2n  Co2n-l  C()2n  i0  A 


The  total  integral  is  then 


I  =  ,(->*  i W.  (-l)n+1,^n"^( — —  (_J _ ]2"->  ♦  - ♦ 

n  n  n  2n-l  InCA^)  (Zn-lKn2""1 


ln(A<)  i  ,  r 

— I  exp[-ie-5)<U)  . 
?o2n  «o2n  Cc  A 


(123) 


By  changing  the  integration  variable  back  to  <,  the  integral  on  the 
right-hand  side  of  (123)  becomes 


/  exp [-  -  e"A]dC  =  j  e — —  d<  =  E^^z)  , 
c  A  k 

m 


(124) 
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where  E \ ( <mz )  is  the  exponential  integral  of  order  1.  For  z<<l/<m,  we 
can  apply  the  small  argument  approximation16  of  E^^z). 


E i  ( x^z )  =  -ln(r<mz). 


Hence,  the  integral  I n  can  approximated  by 


(125) 


T  '  ,  lNn+l  2n-2r  1 

I  =  -1  IT  { - .  I - 

n  2n-l  Lln(A<m) 


1  ,2n-l  2n  r  1 


2n-l  In (A/z ) 


+  In [A/( Tz )  ]  j 
1 n2n  (A/z ) 


,2n-l 


(126) 


for  z.<<z<<1/k  . 

d  m 


From  the  definition  of  A  given  in  (106)  and  z^  given  in  (97), 


we  have 


2  2 

|A|  =!4-z 

4b  < 


(127) 


Since  r  /4=0.8  and  a<b,  the  condition  that  z>>z^  guarantees  that  JA|/z<<1 
and  J 1 n (A/z )  j  >1 .  hence,  it  is  useful  to  note  that  (126)  shows  that 


1  ’2n-l 


I  =  Oli - - —  j 

n  1 n  (A/z  ) 


(128) 


for  z<l/<m. 
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This  means  that  the  series  expansion  of  I(z)  given  in  (116)  is 
dominated  by  I j,  .  From  (126),  we  have 


_l _  _  _2 _  +  1n[A/(rz)] 

1 n ( Ac  m )  ln(A/z)  1n2 (A/z) 


(I2y) 


By  combining  the  last  two  terms,  this  can  be  written 


1 


1 


ln(A<m)  1 n (A/z ) 


n  W/D] 

L  In (A/z ) 


(130) 


The  last  term  in  (130)  is  observed  to  be  the  first  part  of  the  binomial 
expansion 


1 


1 


lr,Aj  1 n (A/z )[ 1+1 n ( 1/r  )/l n (A/z ) ]  ln(A/z) 
rz 


1  ^  _  ln( 1/r ) 


1  n  ( A/z  )J 


+  0 


1 


In  (A/z)' 


-] 


(131) 


Hence,  (130)  becomes 


1 


1  n  (  A«  m  ) 


U  [- 


ln(^— ]  In'  (A/z) 


(132) 


Tz 


Finally,  we  combine  (116),  (128),  and  (132)  to  obtain  the 
expression  for  I(z) 


Hz) 


i 

TnOSJ 


o  [. 


1 


In  (A/z ) 


(133) 
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valid  for  <<z << 1 / <m .  This  result  reduces  to  the  result  in  Section  4  of 
Reference  3  in  the  limit  as  because  ^-*-0  in  that  limit. 


With  the  application  of  (133)  in  (114),  the  portion  Bp(z)  of  the 
branch-cut  integral  can  be  approximated  by 


n  /  \  '  2lTi  ik2z  r  1 
B  (z )  =  — -2—  e  z  1 

P  bk  ik  2 


1 


In  (—)  ln(A<m) 

rz 


] 


The  total  branch-cut  integral  can  then  be  written 


B(Z)  =  eikzZ  [——  -  Q]  , 

bk,k2  1  n  (A — ) 

rz 


where 


(134) 


(135) 


< 

2  m 


„  -KZ  . 
<e  d< 


Q  =  v  +  (kja)  /  . - — - - - - - -  (136) 

1  n  ( A <rn)  0  [Wk2+ak  j<ln(A<)  ]  +  (ak^*) 


for  z,<<z<<1/<  . 
d  m 


The  contribution  of  the  first  term  in  the  brackets  in  (135)  is 
the  asymptotic  value  of  the  branch-cut  integral  for  a  perfect  conductor. 
Hence,  the  contribution  of  the  second  term  is  the  change  introduced  by  the 
fi nite-  conducti vity  of  the  wire.  It  is  nearly  independent  of  z  for  small 
z  and  tends  to  zero  in  the  limit  as  We  determine  its  value 

numerically  for  a  case  of  interest  in  the  next  section. 


SECTION  4 

CURRENT  DISTRIBUTION  ALONG  THE  WIRE 


In  this  section,  we  combine,  summarize,  and  discuss  the  analyti¬ 
cal  results  obtained  in  the  previous  sections  and  present  numerical  examp¬ 
les.  Throughout  the  section,  the  source  is  taken  to  be  on  the  wire  by 
setting  b  equal  to  a. 

4.1  SUMMARY  OF  ANALYTICAL  RESULTS 

The  current  produced  by  a  6  function  voltage  generator  of 
strength  V0  located  on  the  wire  at  z=0  can  be  written  for  positive  z  in 
the  form 


■<z>  ■  W2’ 


+  I  (z ) 

space 


by  combining  (36)  with  (60). 


(137) 


The  first  term  in  (137), 

I  (z)  =  -V 0  kZ2  o^iri^R  (z)  ,  (138) 

mode  j  j 

is  the  contribution  of  the  complex  TM  modes  guided  by  the  wire  with  wave- 
numbers  hj  that  satisfy  the  modal  equation  (43).  All  of  these  modes 
except  one  (known  as  the  principal  mode)  are  damped  out  immediately  away 
from  the  generator.  The  wavenumber  hg  of  the  principal  mode  is  very  close 


to  kz  (the  wave  number  of  the  material  surrounding  the  wire)  and  it  can 

,  2  2  1/2 

lie  on  either  Riemann  sheet  of  the  square  root  ct2  =  (k 2-h o ) •  h0 

lies  on  the  upper  sheet,  the  principal  mode  contributes  to  the  current  in 
the  wire  and  the  mode  is  said  to  have  been  excited.  If  ho  is  on  the  lower 
sheet,  the  mode  does  not  contribute  to  the  current  (Ro  is  then  taken  to  be 
zero)  and  the  mode  is  not  excited. 

When  it  is  excited,  the  current  produced  by  the  principal  mode 
is  given  to  a  good  approximation  for  highly  conducting  wires  by  (76), 

I  (z)  =  _  ,  (139) 

mode  h0k2j21n(I^-)  +  1] 

2i 

2  ~ 

where  a2  =  2k2(k2-h0).  If  we  set 


h0  =  k2+<S  (140) 

in  the  expression  for  a2  and  h0  =  k2  in  the  denominator  of  (121),  then 
(139)  can  be  written 


i  (2)  :  2'’vo^°ietho;  ,  (hi) 

mode  k1[ln(A5)+l] 

where  A  is  given  by  (106). 

In  order  to  evaluate  (141)  for  any  specific  case,  it  is 
necessary  to  solve  the  modal  equation  (43)  numerically  to  determine  h0 
(and  hence  6).  This  is  necessary  also,  in  order  to  determine  whether  or 
not  the  mode  is  excited  by  determining  whether  ho  is  on  the  upper  or  lower 
Riemann  sheet.  It  has  been  found  numerically  that  the  imaginary  part  of 


3> 
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h0  is  always  slightly  larger  than  k^'.  This  means  that  the  attenuation  of 
the  principal  mode  is  slightly  greater  than  that  of  a  plane  wave  in  the 
surrounding  material. 

The  second  term  in  (137), 

*space(z)  =  Vok2ajaB(z)  (142) 

is  the  contribution  of  an  electromagnetic  wave  (known  as  a  space  wave) 
that  radiates  from  the  generator  out  into  the  material  surrounding  the 
wire.  The  wave  number  of  this  wave  is  equal  to  that  of  the  surrounding 
material  but  its  amplitude  decreases  with  distance  from  the  generator 
because  of  the  geometrical  spreading  of  the  wave  as  it  propagates. 

The  distribution  of  current  along  the  wire  is  given  in  (142)  in 
terms  of  the  branch-cut  integral  B(z).  If  z >>z(j=  2 jtc2  | a2  so  that  the 

small  argument  approximations  can  be  used  to  approximate  the  Hankel 
functions  in  the  integrand,  B(z)  is  given  by  (111)  with  b=a.  The 
resulting  expression  for  the  current  is 


I  (z)  =  -2w1V0k2o1a2e1k2Z/" - - -  ,  (143) 

sPace  o  [ Wk2 +ak1»c  1  n (A<  )]  +(ak1iric ) 

with  W  given  by  (10y). 

If  z»zc  =  ( h o -k 2 )  1  ,  then  Ispace(z)  can  be  approximated  by 
(93)  (with  aj  =  kj )  as 

I  (;)  :  4  ?’Vl  [1  ]2l-<k,^  eik2*  ■  (144) 

space  k2  W  z 


This  is  the  current  that  would  be  generated  in  the  far  field  of  an  elec¬ 
tric  dipole  located  at  the  origin  and  radiating  into  a  homogeneous 
material  with  wave  number  k2.  The  reason  that  the  current  falls  off 
faster  than  1/z  is  that  it  is  driven  by  the  longitudinal  electric  field 
rather  than  the  transverse  electric  field  of  the  dipole. 

Since  the  distance  zc  is  very  large  for  a  highly  conducting 
wire,  it  is  desirable  for  applications  to  study  the  properties  of  the 
current  much  closer  to  the  source.  In  the  region  zc|<<z<<l/icrn,  the 
space-wave  current  can  be  written  [by  combining  (142)  and  (135)]  as 


I  (z)  =  2TriVoai  eik2Z[ 


1 


space 


ln[A/(rz)] 


-  Q]  , 


(145) 


where  Q  is  given  by  (136).  The  positive  number  <m  is  the  smallest  value 
of  k  for  which  the  term  Wk2  in  the  denominator  in  (143)  can  be  neglected. 
The  first  term  on  the  right-hand  side  of  (145)  gives  the  current  that 
would  result  if  the  wire  were  a  perfect  conductor.  The  second  term  gives 
the  change  introduced  by  the  finite  conductivity  of  the  wire.  It  tends  to 
zero  as  .  Moreover,  Q  is  nearly  independent  of  z  for  z<<l/<m. 

Hence,  Ispace^2)  differs  from  the  perfect-conductor  current  distribution 
by  having  a  component  that  does  not  fall  off  due  to  geometrical  spreading 
of  the  space  wave  in  this  region. 

The  primary  difference  between  the  current  distributions  due  to 
the  principal  mode  and  the  space  wave  is  that  the  variation  of  the  modal 
current  with  z  is  determined  by  only  its  wavenumber  ho  whereas  the  space- 
wave  current  varies  with  wavenumber  k2  with  an  additional  fall  off  in 
amplitude  due  to  geometrical  spreading  of  the  space  wave.  Since  ho  has 
slightly  larger  imaginary  part  than  k2,  the  space-wave  current  will  always 
dominate  the  modal  current  for  very  large  z.  Closer  to  the  source. 


/* 
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however,  both  currents  can  be  important.  Since  ho  is  very  nearly  equal  to 
k2  and  since  the  geometrical  fall-off  of  the  space  wave  is  very  slow  in 
the  near  region,  the  variation  of  the  currents  with  z  can  be  nearly  the 
same.  The  ratio  of  the  amplitudes  of  the  two  current  is  approximated  by 

=  [ In (A<$)+1  ]  U  -  Q]  e-6"*  (146) 

Imode  ln[A/(Tz)] 

for  Zd«z«l/<m. 

4.2  NUMERICAL  RESULTS 

In  order  to  obtain  more  details  concerning  the  current 
distribution  along  the  wire  and  to  check  the  approximations,  it  is 
necessary  to  evaluate  Imode(z)  and  *space(z)  numerically.  The  first 
is  computed  from  (141)  once  h0  has  been  determined  by  solving  numerically 
an  approximate  form  of  the  model  equation  that  is  valid  for  the  principal 
mode.  When  we  apply  the  approximation  =  k i  and  the  small  argument 
approximations  of  the  Hankel  functions  in  (42),  the  modal  equation  (43) 
reduces  to 

Wk2  +  i  a2ak 1 1 n (  — =  0  .  (147) 

This  reduces  to  the  equation  treated  by  Sommerfeld  to  study  the  modes  of  a 
wire  if  W  is  taken  to  be  one.  Unfortunately,  the  numerical  procedure 
applied  by  Sommerfeld  to  solve  (147)  does  not  converge  for  parameter 
values  of  interest  in  EMP  problems  when  the  correct  value  of  W  is 
applied.  Consequently,  we  applied  Muller's  method  17  to  solve  (147) 
numerically  for  a2.  Then  6  was  evaluated  by 


according  to  (14U)  and  the  equation  just  above  it. 


ISpace  (2)  was  computed  numerically  from  (143)  by  using  the 
trapezoidal  rule  to  evaluate  the  integral.  The  approximations  of 
I$pace(z)  given  in  (144)  and  (14b)  also  were  evaluated  numerically  and 
compared  with  the  value  obtained  from  (143).  In  order  to  apply  (145),  Q 
given  by  (136)  was  computed  using  the  trapezoidal  rule.  was  taken  to 
be  the  value  of  k  at  which  the  error  in  approximati ng  the  integrand  in 
(143)  (by  neglecting  the  term  Wk2 )  is  0.1%  of  the  true  integrand. 

The  values  of  the  parameters  used  for  the  computations  are 
presented  in  Table  1.  The  wire  parameters  correspond  to  a  copper  wire 
with  resistance  of  R  =  10"3  ohm/m.  Additional  computations  were  made  for 
a=7 .44U31 19  x  10“3m  and  a  =  7.44U3119  x  10-4m  which  correspond  to  R  =  10-4 
ohm/m  and  R  =  10-2  ohm/m  respectively.  Since  the  wire  conductivity  is  so 
high,  the  term  involving  e t  in  the  expressions  (3)  and  (4)  for  kL  was 
neglected. 


Table  1.  Parameters  used  in  the  computations. 

oj  =  5.75  x  107  mho/m 

o2  =  10"  mho/m 

z2  =  2000  eg 

u>  =  10V1 

a  =  2.3528332  x  10-3m 

When  the  modal  equation  was  solved,  it  was  found  that  the  pole 
corresponding  to  the  principal  mode  lies  on  the  lower  Riemann  sheet  for 
all  three  values  of  a.  This  means  that  the  residue  of  the  pole  does  not 
contribute  to  the  current  since  the  sum  in  (138)  is  the  sum  of  the 


residues  of  the  poles  in  the  upper  half  of  the  complex  h  plane  on  the 
upper  Riemann  sheet.  Consequently,  we  have 

I  mode  =  0  (149) 

for  these  cases. 

It  was  found  by  increasing  the  frequency  in  the  R  =  10~3  ohms/m 
case  that  the  pole  crosses  the  branch  cut  onto  the  upper  Riemann  sheet 
when  a)  reaches  2.60  x  104  s_1  .  At  that  frequency,  the  pole  begins  to 
contribute  an  amount  given  by 

Imode(z)  =  "  V0(.O33+.025i)eik2z  (lbU) 

amps  at  z  =  1000m.  As  the  pole  crosses  the  branch  cut,  the  branch-cut 
integral  increases  abruptly  by  minus  the  amount  given  in  (150)  so  that  the 
total  current  does  not  change  discontinuously. 

Computations  also  were  made  for  oj  =  5.75  x  1012  mho/m  to  see  if 
the  results  agreed  with  the  asymptotic  theory  of  a  perfectly  conducting 
wire.  Agreement  was  found  to  at  least  three  significant  figures  for  z 
less  than  about  10,000  m.  Beyond  that  distance,  the  perfect-conductor 
theory  begins  to  break  down  for  this  value  of  oj . 

In  order  to  display  the  numerical  results,  it  is  convenient  to 
factor  out  the  propagation  term  exp(ik2z)  and  the  strength  Vq  of  the 
excitation.  Hence,  Figures  3  and  4  respectively  show  the  real  and 
imaginary  parts  of  Iq  (z  )  defined  by 
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plotted  as  a  function  of  z  on  a  log-log  plot  for  the  parameter  values 
given  in  Table  1  (R  =  10"3  ohm/m). 

The  data  shown  in  these  figures  do  not  agree  with  either  of  the 
approximations  of  ISpace  (z)  presented  in  (144)  and  (145).  To  see  why, 
we  now  determine  numerically  the  regions  of  validity  of  the 
approximations.  Both  approximations  require  that  z>>zc j  in  order  for  the 
smal 1 -argument  approximations  of  the  Hankel  functions  to  be  applicable. 
Evaluation  of  (97)  for  the  R  =  10"3  ohm/m  case  shows  that 
Z(j=1.2  x  10"8  m.  Since  this  value  is  so  small,  we  conclude  that  the 
above  condition  will  be  satisfied  for  any  point  of  interest  along  the 
wi  re. 

The  other  conditions  for  validity  of  the  approximations  are  much 
more  restrictive.  In  particular,  in  order  for  (145)  to  be  valid  with  Q 
defined  by  (136)  we  require  that  z<<l/<m.  Numerical  computation  of  the 
terms  in  the  denominator  of  the  integrand  shows  that  =  0.83  m_1  for 
the  R  =  10-3  ohm/m  case.  (See  the  second  paragraph  above  Table  1  for  the 
definition  of  <m).  This  means  that  we  require  that  z<<1.2  m  in  order 
for  (145)  to  be  useful.  Numerical  evaluation  of  the  integral  expression 
for  Ispace  (z)  ancf  Q  for  very  small  z  shows  that  (145)  agrees  with  the 
true  value  of  ISpace  (z)  to  within  2  percent  for  z<0.1  m.  Hence,  the 
approximation  provides  a  useful  check  of  our  numerical  evaluation  of  the 
integral  but  it  does  not  provide  an  analytic  approximation  of  I$pace  (z) 
in  a  region  of  interest. 

Similarly,  in  order  for  (144)  to  be  useful,  we  must  have 

ii  ll  _  i  . 

z>>zc  =  (ho  -  kz)  .  For  the  R  =  10-  ohm/m  case,  it  is  found  that 
zc  =  7.3  x  103  m.  Numerical  evaluation  of  the  integral  expression  for 
Ispace  (z )  shows  that  the  approximation  given  in  (144)  is  good  to  within 


23  percent  for  z>10,000  m  and  to  within  3.5  percent  for  z>106  m.  Again, 
we  find  that  the  approximation  provides  a  valuable  check  on  the  numerical 
evaluation  of  the  integral  for  Ispace  (z)  but  it  does  not  provide  an 
analytic  expression  that  is  useful  in  the  region  of  interest  along  the 
wi  re. 

The  range  of  validity  of  the  approximations  can  be  better  for 
other  values  of  wire  resistance.  For  R=10“2  ohm/m,  the  value  of  zc  is 
found  to  be  770  m.  Comparison  of  1 0 ( z )  with  the  approximation  correspond¬ 
ing  to  (144)  is  shown  in  Figure  5.  It  is  seen  that  the  approximation  is 
useful  for  z  beyond  about  lOKm. 

For  R=10-4  ohm/m,  is  found  to  be  0.069  m"1.  Hence,  (145) 
is  useful  for  z<<14.5  m.  Figure  6  shows  the  comparison  of  the  real  part 
of  I0(z)  with  the  approximation  corresponding  to  (145).  Agreement  is 
reasonable  for  z  up  to  about  10  m. 

Figure  6  shows  also  the  real  part  of  the  asymptotic  approxima¬ 
tion  of  Io(z)  for  a  perfect  conductor.  The  result  is  rather  surprising. 
Even  at  very  small  distances  from  the  source  with  wire  resistance  as  low 
as  R=10-4  ohm/m,  the  perfect-conductor  result  is  not  very  useful.  The 
reason  for  the  difference  is  easily  seen  by  comparing  the  numerical  values 
of  the  integrands  in  the  expressions  for  Ispace  (z)  for  the  two  cases. 

The  absolute  values  of  the  integrands  are  shown  in  Figure  7  as  a  function 
of  distance  up  the  branch  cut  for  z=100  m.  It  is  seen  that  the  primary 
effect  of  neglecting  the  finite  resistance  of  the  wire  is  to  greatly 
increase  the  magnitude  of  the  integrand  near  the  branch  point.  For  the 
finite  resistance  case,  the  integrand  is  negligible  along  the  portion  of 
the  branch  cut  that  is  below  the  pole.  For  the  perfect  conductor, 
however,  that  portion  of  the  branch  cut  provides  a  significant  contribu¬ 
tion  to  the  integral  which  increases  the  absolute  value  of  the  resulting 
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In  order  to  display  the  effect  that  finite  resistance  has  on  the 
current  distribution  along  the  wire,  let  Ipc(z)  be  the  asymptotic 
approximation  of  Iq(z)  for  a  perfect  conductor  defined  by 


I  (z)  =  - 

PC  k2  ln[A/( Tz ) ] 

and  let  Ir(z)  be  defined  by 

IQ(z)  =  I pc (z)  +  I R ( z ) 


(152) 


(153) 


for  the  resistive  wire*.  With  this  definition  Ir(z)  was  determined 
numerically  for  the  R=10-3  ohm/m  case  by  numerical  evaluation  of  Iq(z) 
and  IpQ ( z ) .  The  results  are  shown  in  Table  2  where  A0,  Apc,  Ar  and 
4>0 »  'J’pc*  4>r  are  respectively  the  amplitudes  in  amp/m  and  phases  in 
radians  of  Io(z),  Ipc ) »  and  Ir(z).  It  is  seen  that  Ir(z)  is  a 
very  slowly  varying  function  of  z  over  a  very  large  range.  Moreover,  it 
is  seen  that  Ir(z)  is  similar  in  magnitude  to  I pq (z )  but  nearly 
opposite  in  phase.  As  a  result,  the  amplitude  of  the  current  is  reduced 
from  the  perfect-conductor  case  by  a  substantial  amount  which  is  given  the 
last  column  in  Table  2.  The  results  show  that  the  decrease  in  current 
amplitude  is  nearly  independent  of  distance  from  the  source  over  the 
entire  range  of  interest.  Since  the  current  in  a  perfect  conductor 
decreases  with  distance  from  the  source,  the  relative  change  in  current 
due  to  wire  resistance  increases  with  distance. 


*  Note  that  for  z d<<z<<  1/Km,  1r(z)  is  approximately  proportional 
to  Q  according  to  (145). 


66 


SECTION  5 
CONCLUSIONS 


The  results  of  this  study  have  exposed  the  reason  that 
transmission-line  theory  does  not  describe  accurately  the  behavior  of 
currents  generated  by  localized  sources  on  uninsulated  wires.  The  current 
is  dominated  by  an  electromagnetic  space  wave  that  radiates  energy  out 
into  space.  Since  transmission-line  theory  applies  to  IEM  modes  guided  by 
the  transmission  line,  it  does  not  describe  the  behavior  of  space  waves. 
Although  the  wire  does  support  a  guided  TM  mode  that  is  approximately 
transverse  electric  and  hence,  could  be  described  by  transmi ssi on-1 i ne 
theory,  the  frequencies  of  interest  are  too  low  for  the  mode  to  be 
excited.  For  frequencies  above  several  KHz,  the  mode  begins  to  be  excited 
but  it  still  makes  only  a  minor  contribution  to  the  total  current. 

A  simple  analytic  approximation  is  available  for  the  current 
distribution  along  the  wire  in  the  region  of  interest  provided  that  the 
wire  is  a  perfect  conductor.  This  study  has  provided  the  numerical  values 
of  the  correction  term  that  must  be  added  to  account  for  the  finite 
resistance  of  the  wire  for  typical  conductivities.  It  was  found  that  the 
correction  term  is  nearly  constant  along  the  wire  and  that  it  is  similar 
in  magnitude  but  nearly  opposite  in  phase  from  the  perfect-conductor 
result.  Consequently,  an  important  effect  of  finite  wire  resistance  is  to 
reduce  the  current  amplitude  by  a  substantial  amount  that  is  nearly 
independent  of  distance  from  the  source.  The  relative  difference  in  the 

q 

current,  however,  increases  with  distance  from  the  source.  For  R  =  10- 
ohm/m,  the  current  amplitude  is  about  20  percent  of  the  perfect-conductor 
result  for  z=lm  and  about  5  percent  for  z=1000m. 
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AVCO  Systems  Division 
ATTN:  Library 


Franklin  Institute 

ATTN:  R.  Thompson 
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General  Dynamics  Corp 

IRT  Corp 

ATTN 

Research  library 

ATTN 

:  B.  Williams 

ATTN 

:  N.  Rudie 

General  Dynamics  Corp 

ATTN 

Research  Library 

ITT  Corp 

ATTN 

:  A.  Richardson 

General  Electric  Co 

ATTN 

:  Technical  Library 

ATTN 

J.  Andrews 

JAYCOR 

General  Electric  Co 

ATTN 

W.  Radasky 

ATTN 

C.  Hewison 

ATTN 

D.  Higgins 

General  Electric  Co 

JAYCOR 

ATTN 

Technical  Library 

ATTN 

R.  Poll 

ATTN 

R.  Schaefer 

General  Research  Corp 

3  cy  ATTN 

Technical  Info  Office 

JAYCOR 

ATTN 

Library 

Georgia  Institute  of  Technology 

ATTN 

Res  &  Sec  Coord  for  H.  Denny 

JAYCOR 

ATTN 

E.  Wenaas 

Grumman  Aerospace  Corp 

ATTN 

R.  Stahl 

ATTN 

L-01  35 

Johns  Hopkins  University 

GTE  Communications  Products  Corp 

ATTN 

P.  Partridge 

ATTN 

J.  McElroy 

Kaman  Sciences  Corp 

GTE  Coninunications  Products  Corp 

ATTN 

A.  Bridges 

ATTN 

H.  Gel  man 

ATTN 

F.  Shelton 

ATTN 

J.  Waldron 

ATTN 

W.  Rich 

ATTN 

A.  Murphy 

ATTN 

N.  Beauchamp 

ATTN 

D.  Flood 

ATTN 

M.  Snow 

Kaman  Sciences  Corp 

ATTN 

F.  Harmond 

ATTN 

E.  Conrad 

GTE  Government  Systems  Corp 

Kaman  Tempo 

ATTN 

L.  Lesinski 

ATTN 

DAS  I AC 

ATTN 

R.  Steinhoff 

Kaman  Tempo 

Harris  Corp 

ATTN 

DAS  I  AC 

ATTN 

A.  Strain 

ATTN 

W.  Hobbs 

ATTN 

V  Pres  &  Mgr  Prgms  Div 

ATTN 

R.  Rutherford 

ATTN 

W.  McNamara 

Hazel  tine  Corp 

ATTN 

J.  Okrent 

Litton  Systems,  Inc 

ATTN 

MS  64-61,  E.  Eustis 

Honeywell,  Inc 

ATTN 

R.  Johnson 

Litton  Systems,  Inc 

ATTN 

S&RC  Library 

ATTN 

J.  Skaggs 

Honeywell,  Inc 

Litton  Systems,  Inc 

ATTN 

S.  Grafe 

ATTN 

J.  Moyer 

ATTN 

W.  Stewart 

Lockheed  Missiles  &  Space  Co,  Inc 

Hughes  Aircraft  Co 

ATTN 

B.  Kimura 

ATTN 

A.  Narevsky 

ATTN 

L.  Rossi 

ATTN 

S,  Taimuty 

Hughes  Aircraft  Co 

ATTN 

H.  Thayn 

ATTN 

K.  Walker 

ATTN 

D.  Nishida 

ATTN 

CTDC  6/E  110 

Lockheed  Missiles  &  Space  Co,  Inc 

1 1 T  Research  Institute 

ATTN 

Technical  Info  Center 

ATTN 

I.  Mindel 

ATTN 

J.  Bridges 

Lutech,  Inc 

ATTN 

F.  Tesche 

Institute  for  Defense  Analyses 

ATTN 

Tech  Info  Services 

Maqnavox  Govt  ft  Indus  Electronics  Co 

ATTN 

W.  Richoson 
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Martin  Marietta  Corp 

ATTN:  M.  Griffith 
ATTN:  J.  Casalese 

McDonnell  Douglas  Corp 
ATTN:  M.  Potter 
ATTN:  W.  McCloud 

McDonnell  Douglas  Corp 

ATTN:  S.  Schneider 

ATTN:  Tech  Library  Services 

McDonnell  Douglas  Corp 


ATTN 

T. 

Ender 

Metatech  Corp 

ATTN 

w. 

Radasky 

Mission  Research 

Corp 

ATTN 

C. 

Longmire 

ATTN 

EMP  Group 

2  cy  ATTN 

G. 

Sherman 

5  cy  ATTN 

Document  Control 

Mission  Research 

Corp 

ATTN 

M. 

Scales 

ATTN 

A. 

Chodorow 

Mission  Research 

Corp 

ATTN 

W. 

Stark 

ATTN 

J. 

Lubel 1 

ATTN 

w. 

Ware 

Mitre  Corp 

ATTN:  M. 

Fitzgerald 

Norden  Systems,  Inc 

ATTN:  Technical  Library 
ATTN:  D.  Longo 

Northrop  Corp 

ATTN:  Rad  Effects  Grp 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 

Palisades  Inst  for  Rsch  Services,  Inc 
ATTN:  Records  Supervisor 

Physics  International  Co 

ATTN:  Document  Control 

R&D  Associates 

ATTN:  Document  Control 
ATTN:  C.  Knowles 
ATTN:  W.  Karzas 
ATTN:  C.  Mo 
ATTN:  P.  Haas 

Rand  Corp 

ATTN:  Lib-D 
ATTN:  P.  Davis 

Rand  Corp 

ATTN:  B.  Bennett 

Raytheon  Co 
ATTN 
ATTN 
ATTN 
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RCA  Corp 


ATTN: 

D.  O'Connor 

ATTN: 

L.  Minich 

RCA  Corp 

ATTN: 

G.  Brucker 

Rockwell  International  Corp 

ATTN: 

N.  Rudie 

ATTN: 

J.  Erb 

ATTN: 

V.  Michel 

ATTN: 

D/243-068,  031-CA31 

ATTN: 

G.  Morgan 

Rockwell  International  Corp 

ATTN:  8-3  Div  Tic,  BAOB 


Rockwell  International  Corp 
ATTN:  F.  Shaw 

Rockwell  International  Corp 
ATTN:  B.  White 

Sanders  Associates,  Inc 

ATTN:  R.  Despathy 

Science  &  Engry  Associates,  Inc 
ATTN:  V.  Jones 

Science  Applications  Inti  Coro 
ATTN:  E.  Parkinson 

Science  Applications,  Inc 
ATTN:  E.  O’Donnell 

Science  Applications,  Inc 
ATTN:  W.  Chadsey 

Sidney  Frankel  &  Associates 
ATTN:  S.  Frankel 

Singer  Co 

ATTN:  Technical  Info  Center 
Sperry  Corp 

ATTN:  Technical  Library 

Sperry  Corp 

ATTN:  M.  Cort 

Sherry  Flight  Systems 
ATTN:  D.  Schow 

SRI  International 

ATTN:  A.  Padgett 

SRI  International 

ATTN:  A.  Whitson 
ATTN:  M.  Tarrasch 

Teledyne  Brown  Engineering 
ATTN:  F.  Leopard 
ATTN:  J.  Whitt 

Texas  Instruments,  Inc 
ATTN:  D.  Manus 
ATTN:  Technical  Library 


H.  Flescher 
M.  Nucefora 
B.  Schupp 


Transients  Limited  Coro 
ATTN:  D.  Clark 


